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We investigate the coherence properties of thermal atoms confined in optical dipole traps where 
the underlying classical dynamics is chaotic. A perturbative expression derived for the coherence of 
the echo scheme of [Andersen et. al., Phys. Rev. Lett. 90, 023001 (2003)] shows it is a function of 
the survival probability or fidelity of eigenstates of the motion of the atoms in the trap. The echo 
coherence and the survival probability display "system specific" features, even when the underlying 
classical dynamics is chaotic. In particular, partial revivals in the echo signal and the survival 
£f) ■ probability are found for a small shift of the potential. Next, a "semi-classical" expression for the 

averaged echo signal is presented and used to calculate the echo signal for atoms in a light sheet 
wedge billiard. Revivals in the echo coherence are found in this system, indicating they may be a 
generic feature of dipole traps. 
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■ I. INTRODUCTION 

' The capability to control and manipulate quantum systems that has been developed in recent years holds promise 
Q\ , for applying quantum systems to information processing [1]. This possible application has drawn much attention 
to the subject of decoherence/dephasing, where the concept of "fidelity" and "survival probability" plays a major 
role. Some of the possible quantum information schemes involve superposition states of the internal states of trapped 
neutral atoms or ions. Interaction with the trapping potential often causes these superposition states to dephase, but, 
in analogy with the spin echo [2], a coherence "echo" can be used to preserve the coherence [3] [4] [5]. Under certain 
conditions the echo coherence of trapped atoms is limited by the quantum dynamics of the atoms in the trap, thereby 
enabling study of quantum dynamics by echo techniques, but also stressing that the quantum dynamics of the atoms 
in the trap is of great importance to the coherence of their internal states [4] . 

In this paper the dephasing of trapped atoms in coherent superposition states of their internal states is investigated. 
A special focus is put on neutral atoms trapped in optical dipole traps, where the underlying classical dynamics is 
chaotic namely atom-optics billiards [6] [7]. This is done partly because chaotic billiards where used in the past as 
a simplified model to describe decoherence [8], and partly because they can have very long coherence times [9]. The 
atoms are considered " hot" , meaning that their temperature is much larger than the mean level spacing of the trap, 
so they occupy high excited states in the trap, and an output of the experiment is an ensemble average over many 
occupied states. 

The strength of the dipole potential depends on the internal state of the atom. This causes atoms initially in coherent 
superposition states of their internal degree of freedom but with different external initial conditions to acquire different 
phases thereby loosing the macroscopic coherence. This dephasing can be suppressed by echo techniques, where each 
part of the superposition state spends an equal amount of time in each potential, thereby ensuring that the phases 
acquired are equal [4] [5]. However, since the external potentials arc slightly different (we denote this difference "the 
perturbation"), also the dynamics depends on the internal state, and loss of coherence from this is not cancelled by 
echo techniques. For far-off -resonance optical dipole traps such a perturbation is small, so a perturbative result is 
of great interest. In section III a perturbative expression for the echo signal is derived, and confirmed numerically, 
showing that the echo signal can be written as a function of the survival probability or fidelity of eigenstates of the 
motion of the atoms in the trap. Based on this perturbative result the long-time echo coherence of atomic ensembles 
is predicted. 

In the perturbative limit the echo coherence and the survival probability display "system specific" features, despite 
the fact that the underlying classical dynamics is ergodic [10]. This is demonstrated in section IV by showing partial 
revivals in the echo signal and the survival probability for a perturbation consisting of a small spacial shift of the 
potential. In section V a " semiclassical" perturbative expression for the averaged echo coherence is derived. Two 
billiard systems are considered throughout the paper. First, an "idealized" system - a shift perturbation of the 
hard-wall Bunimovich billiard [12] where eigenstates can be computed efficiently [13] is used to confirm the validity 
of the perturbative and semiclassical approximations by comparing them to a full quantum calculation (sections IV 
and V). These approximations are then used in section VI to calculate the echo coherence for micro wave (MW) 
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spectroscopy of atoms in an experimental realizable atom optics billiard, namely the light-sheet wedge billiard [6] [14], 
where revivals in analogy with those in the hard-wall stadium are found. This suggest that revivals of coherence may 
be a generic feature of optical dipole traps. 

II. MICRO WAVE SPECTROSCOPY AND QUANTUM DYNAMICS OF TRAPPED ATOMS 

In this section pulsed micro wave spectroscopy on optical trapped alkali atoms is considered. 

The ground state of alkali atoms contains two hyperfine states with an energy splitting of Ehf — ^hf, each 
containing several magnetic substates. The hyperfine energy splitting in Caesium atoms is the foundation of modern 
time scales, so MW spectroscopy of this transition in an isolated environment is of great importance and is therefore 
a very advanced field. However, many quantum information processing schemes require trapped atoms or ions, and 
if neutral alkali atoms are placed inside an optical trapping potential, then the energy levels are changed. Since the 
optical dipole potential is inversely proportional to the detuning of the trap laser from resonance [15] and the detuning 
will differ by the hyperfine splitting, then there is a slight difference in potential for atoms in different hyperfine states, 
even if the states considered have the same matrix element for the dipole interaction with the trapping light. For such 
two states (denoted |f) and | J.)) the Hamiltonian considered can be written as, 

H = H l ||> (|| + (iff + Ehf) It) (T| = (f^ + V i (x)) ||) (|| + (£- + V T (x) + E HF ^J |T) (Tl , (1) 

where p is the atomic center of mass momentum and Vf [Vf] the external potential for an atom in state ||) [If)]. 
Typically V is much smaller than Ehf and (Vf — VjJ /Vf <~ 1CP 3 — 10 -5 [4] [5]. The eigenenergies of this Hamiltonian 
consist of two manifolds (belonging to ||) and |f )) separated in energy by Ehf (see Fig. 1). 

The initial step in MW spectroscopy is to prepare all the atoms in one of the spin states involved (||)), so their total 
wave function can be written as | J.) \ipi) = |l, ipi), where \ipi) represents the motional (external degree of freedom) 
part of their wave function. If atoms prepared in | j,, ipi) are irradiated with a MW field close to resonance with ojhf, 
then transitions between eigenstates belonging to different manifolds can be driven. The matrix elements for these 
transitions is given by C n k = (&f | nj.) x Mf_>f where Mf_>f is the free space matrix element for the internal state 
transition, and (&f | n{) is the overlap between \n{), the initial eigenstate of iff, and \k-\), the final eigenstate of iff 
(here the very small momentum of the MW photon is neglected). Since in general each of the eigenstates in the 
|T)-manifold is coupled to many eigenstates in the |T)-manifold by the MW-field, there is no simple prediction of the 
wave function after MW irradiation. However, in the simple limit of very short MW-pulses the external wave function 
is left unchanged and only the internal state of the atom is changed. Such a short pulse that puts atoms initially in 
\i,ipi), into the coherent superposition state \l,ipi) + ITiV'i) is called a ^-pulse and a pulse that flips the spin 
to i\1,ipi) is called a 7r-pulse. 

In the widely used Ramsey spectroscopy technique [16], atoms initially prepared in ||,V'i) (assumed to be a pure 
state) and irradiated with a ^-pulse to generate the wave function -4^ \l,ipi) + -75 |T>^i)- After some time r this 

state will, in the rotating frame of the MW field, evolve into -4^ exp (^—i^rj ||, ipi) + ^ cxp ^— i (j^ + Aj rj |T, ipi), 

where A = uhf — ^mw is the MW-detuning. Then the atoms are irradiated with a second |-pulse generating the 

wave function \ (exp (-«-^r) - cxp (-i + A^j ||, ipi) + \ (^cxp (-i + + cxp (-^ir 7 ")) IT.V'i)- 

The population of state \\) is then, 

cos(Ar + <A i (T))^ , (2) 

where (pi (r) is a phase depending on the initial state \ipi) and the dynamics of it in the trap. By choosing A 
sufficiently large, compared to variations in (pi (r) Eq. 2 yields the well known cosinusodial Ramsey-fringes as a 
function of At, but with a fringe contrast reduced to the square root of the fidelity of the external motion in the trap 

2 

. In the case where \ipi) is an eigenstate of iff, the fidelity can be written 




= Psrv and is therefore often denoted the survival probability of \ipi) in iff. Because 
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schemes where it measures the overlap between the desired state e *( h ' \ipi) and the actual state e l ( ^ ) T and 
because it can be used as an indicator of quantum chaos, its decay have been an immensely studied topic throughout 

However, experiments are often not performed with atoms initially prepared in a pure initial state but with 
atoms in a thermal mixture of states. Since fa (r) depends on the initial state, then the fringe contrast in a Ramsey 
experiment will decay rapidly, not due to a decay in the fidelity of the individual states thermally populated, but 
because the cosine terms from different populated states get out of phase [4] . This rapid loss of coherence is a problem 
in quantum information schemes [5] and prevents the use of Ramsey spectroscopy for studying quantum dynamics [4] . 
The loss in coherence may be overcome by adding an additional potential making = , and thereby decoupling 
the external dynamics from the spectroscopic evolution [18]. This is however not always possible and does not enable 
the use of MW-spectroscopy for dynamical studies. The echo pulse sequence demonstrated in [3] [4] [5] provides a 
more general solution to suppress the inhomogeneous dephasing of superposition states. 

The echo pulse sequence consists of three short pulses (| — n — j), separated by two dark periods of equal duration 
t, after which the population of ||) is measured. In analogy to Eq. 2 the population of state ||) can be written as: 




Fech = 1 indicates complete preservation of coherence and yields Pj = 0, whereas Fech = 0, yielding P] = j 
indicates complete loss of coherence [19]. From Eq. 3 it is seen that the echo amplitude is insensitive to the addition 
of constants to iff and iff, so in contrary to the Ramsey signal (Eq. 2) Pj of Eq. 3 does not depend on u>mw 
and Ehf, but only on the quantum dynamics of the atoms in the potentials. The rapid loss of coherence when a 
thermal mixtures of states is considered has vanished, making echo techniques of interest in experiments on quantum 
information processing [5]. This also means that even for thermal mixtures of states echo spectroscopy can yield 
information on the quantum dynamics of the trapped atoms, as was demonstrated in [4] for more than 10 6 thermally 
populated states. 

III. THE CONNECTION BETWEEN THE ECHO AND THE SURVIVAL AMPLITUDES AND LONG 

TIME BEHAVIOR. 

In this section a second order perturbative expression for the echo amplitude of an eigenstate is derived, showing 
that it can be written as a function of the fidelity or survival amplitude, thereby also predicting the time average of 
the echo amplitude. 

Since thermal atoms can be considered as a mixture of eigenstates, the attention is now turned to the echo amplitude 
of an eigenstate of iff. The echo amplitude of an eigenstate \m\) of iff can be rewritten: 

Fech (H) ,t) = ^ cxp (i (uj l m + w\ - uj l k - wj) r) (m x | Z T ) (Z T | fcf ) (fy | n T ) (n T | m x ) (4) 

l,k,n 

Since adding constants to iff and iff does not change the echo we can without loss of generality set ujj n — = 0. 
Now for atoms trapped in far off resonance optical dipole traps made of linearly polarized light the difference between 
iff and iff is very small, so a perturbative result is of great interest. For a small perturbation (n-f | nj,) ~ 1 — and 
(fcf | ni) is of the order e or less. Rewriting Eq. 4 and keeping only terms up to order e 2 yields: 

Fech {\mi) ,t) ~2 ^ ex P (-"4 T ) l( n T I m l)| 2 + 

exp (m4t) |(n T | mi)\ 2 - 
^exp (i (wj -u^) T ) I <nt \ mi)\ 2 - 
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exp (-« ( W T + w l) T ) |(n T | m 1 >| 2 + 

|(m T | m^l 4 (5) 

This rewriting does not look like a simplification of Eq. 4, but first it is noted that the echo is no longer a function 
of matrix elements like (nj | mi) but only a of their absolute value square, often referred to as the local density of 
states (LDOS) [10]. Second, note that since = uij^ = 0, then typically uj^ — uj^ << ut^, so except for extreme 
long times the approximation uj\ ~ = uj n can be used. For high excited states X) n ^m ex P (* T ) ( n T I m l)| 2 = 

Sn^m ex P ( — * ( w n) T ) I ( n t I m l)| 2 i s a l so a good approximation, for times much shorter than 1/A m ; s with A m ; s the 
mean level spacing. Using these approximations Eq. 5 can be rewritten as: 

Fech ,t) ~ 4 exp(-iw„T) |(n T | m^ 2 - 

^2 exp(-2z^ n r) |(n T | mj.)| 2 + 

IKIm,)! 6 , (6) 

where \(m-\ \ m{)f ~ 1 — 6e^. The quantity ^m^ exp ^— i — rj m^j is denoted the survival amplitude 

(Fsrv), since its length square is the survival probability. Writing the survival amplitude in the same way as Eq. 6 
shows that Eq. 6 yields a simple relation between the echo amplitude and the survival amplitude of an eigenstate in 
the perturbative limit: 

%h(N,t) ~4F S Rv(\m l ),T)-F S Rv(\m l ),2T)-2 (7) 

Eq. 7 indicates that the time average of the echo amplitude is always smaller than the time average of the survival 
amplitude. Hence in the low temperature regime, where the decay of coherence is set by the decay of the survival 
amplitude, and not by the thermal spread of initial conditions, a Ramsey type experiment will yield better coherence 
than an echo experiment (see Eq. 2). This is in agreement with the experimental results published in [3]. 

From Eq. 6 it is seen that for times r where the approximation uj^ ~ uj^ is valid, the constant contribution to the 
echo signal is not given by the infinite long time limit |(mj | mj,)| but by |(to| | m^)| . This is a rather peculiar result 
since |(m| | tij.)| 6 < |(»nj | r/x^ ) | 4 hence it implies that waiting for dephasing to occur on time scales of 1/ (oj^ — w^) 
will increase the ensemble average of the echo coherence. For perturbations that conserve the cigcnspcctrum the 
|(m| | mi)\ 6 echo level persists for all times. 

In [4] we compared the measured long time echo signal for atoms confined in a Gaussian laser beam to a calculation 
based on the ensemble average of the infinite time level of | (to| | } | 4 . However, approximating the Gaussian trap of 
[4] by a harmonic oscillator yields that typically 2n/ (co^ — oj^j ~ 5s, whereas the long time echo level was measured 
when the initial oscillations had died out, typically after r ~ 10ms. Hence, from the perturbative treatment we 
expect to observe the ensemble averaged |(m-f | m^f echo level rather than the \(m^ \ m^f level. We repeated the 
echo- level calculation of [4] using now the assumption that the long time echo signal is given by |(to| | m[}\ , and 
found better agreement with the experimental data published in [4] as compared to the \ (m^ \ ) | 4 calculation, but 
due to the uncertainty of the data points no definitive conclusion can be made [11]. 

Finally, the simplified model presented here predicts that if \(m^ \ m^)f ~ 1 for all occupied states, then perfect 
echo coherence will persist for infinite times. In real systems, there are effects, such as photon scattering or temporal 
fluctuations in the trap depth not taken into account in our model (Eq. 1) causing echo coherence to be lost on long 
time scales [4] [5]. 



IV. EQUATION 7, AND DECAY AND REVIVALS OF ECHO AND SURVIVAL PROBABILITY. 

In this section the validity of Eq. 7 and the behavior of the quantities involved in it, are investigated by numerical 
calculations on a chaotic hard wall Burnimovich stadium [12]. This system is used as a test system, since the 
Burnimovich stadium has become almost a synonym with classical chaotic billiard dynamics, and since the scaling 
method of [13] [20] enables efficient computation of high excited eigenstates of this system. It is demonstrated that 
the echo amplitude and the survival probability show partial revivals at certain times for a spacial shift perturbation. 
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The stadium considered consists of two semicircles of radius 1, connected by two lines of length 2 (see fig. 2). 
We consider the perturbation caused by shifting V\ (x) a small amount along the long axis of the stadium compared 
to Vi (x). This perturbation conserves the eigenspectrum so u\ = uj^. The shift also has direct relevance to the 
"dephasing free" atom optics billiard proposed in [9], and illustrates some of the points of later sections. 

The computation of the echo and the survival probability is done in the straight forward way. A large cluster 
containing <~ 1300 eigenstates with energies around 11000 is computed using the scaling method of [13]. The units h 
=2m=l are used so the eigenstates in the interior of the billiard are solutions to V 2 \ip) = —k 2 and are vanishing 
on the boundary. The transformation matrix T taking states from a representation in eigenstates of to a 
representation in eigenstates of H± is then computed by computing the overlap integral of the eigenstates of 
with those of H± (the length square of the elements of T being the LDOS [10]). The computation of the overlaps is 
efficiently done using the methods of [20] . Since the time evolution of a state in the representation of eigenstates of 
the Hamiltonian under which it evolves is trivial, it is a trivial task to compute the echo and the survival amplitude 
when T is known using Eq. 7. The issue of T not being completely unitary for a hard wall billiard is addressed in 
[10]. 

The structure of T can depend on the perturbation even when the underlying classical dynamics is chaotic [10] [20] 
[21]. In Fig. 3 a plot of the square of the elements of T is shown for the perturbation of a shift of 0.004, where two 
distinct sidebands are seen. This structure induces special behavior of the echo and the survival probability, since an 
eigenstate of H± will show a (partial) revival when evolved in at the time the sidebands have acquired a phase of 
2n. This is demonstrated in Fig 4, where the echo and survival probability for several different eigenstates are shown. 
They all show revivals at the time tr corresponding to the sidebands acquiring a phase of 27r. At longer times the 
revivals from different states get out of phase due to the finite width of the sidebands in Fig. 3. We note that the 
revivals seen in Fig. 4 are not the trivial revivals originating from the system being discrete, but they appears since 
the nature of the difference between H± and induces " selection rules" that determines to which states in iJj an 
eigenstate of is coupled. 

Shown in Fig. 5 is the perturbative approximation Eq. 7 to the echo, together with the exact calculations for a 
shift of 0.001 (Fig. 5a) and 0.008 (Fig. 5b), showing good agreement, even for large perturbations. In 5c the RMS 
deviation between the echo and the approximation (Eq. 7) for the time interval 0-0.1 is shown as a function of the 
shift. It is seen that the deviation grows as e 3 , as expected from Eq. 7 being a second order approximation to the 
echo. This indicates that the other approximations made are valid. 

In Fig. 6 the difference between the exact echo and the approximation (Eq. 7) is shown, for long times, for one of 
the states of Fig. 4. No divergence as a function of time is observed. This is subscribed to the fact that the terms 
neglected in order to obtain Eq. 7 are oscillatory in time. 



V. SEMI-CLASSICAL ESTIMATION OF THE ECHO AMPLITUDE 



In this section the echo amplitude of an ensemble of high excited eigenstates is calculated using classical estimations 
of the LDOS [10]. The method is confirmed by comparison to a full quantum calculation on a hard-wall Burnimovich 
stadium. It is then used in section VI to calculate the echo in a chaotic atom optics billiard, that have been 
demonstrated experimentally in the past [6] [14]. Revivals in analogy with those of the previous section are found to 
appear. 

6 9"'^ 

From the fact that the perturbative limit is considered so |(mj | m\)\ ~ |(mf \m\)\ (where x stands for the 
average of x) it is apparent that Eq.s 5-7 also hold for averages. This means that classical methods [10] [20] [21] can 
be used to estimate the average LDOS and from this the averaged echo amplitude can be computed. 

A detailed description of the computation of the classical estimations of the LDOS is given in [10] [20] [21], so we 
only give a brief review of the steps involved. Writing = Hi+ eH p then first order perturbation theory gives that: 

("T I m l)| 2 — ir^2 — \( n l\ H p \ m l)\ 2 (8) 

For high excited states the band profile of the matrix |(nj H p |rn|)| 2 can be estimated classically as: 

\( ni \H p \ mi )f(AE n , m ) * (^p) , (9) 

where the average is taken over few adjacent states, A m ; s is the mean level spacing and C (u) is determined as 
following. A long classical trajectory x (t) is computed, and — H p is evaluated along it. The time average F m is 
subtracted to yield the fluctuating quantity: 
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F(t) = -H p (x(t))-F m (10) 

The autocorrelation of F (t) is called C (t), and C (co), its Fourier transform, can be computed as the power spectrum 
of F (t). C (lu) is an ensemble average of all classical trajectories of energy E m , but since a chaotic system is considered, 
this can be taken as one long trajectory. The deeper reasons for the validity of Eq. 9, and its good accuracy, have 
been studied in detail [10] [20] [21]. Combining Eq. 9 and 8 yields the expression: 



|(mj | m^)\ 2 is then found as: 

|(m T | mi )\ 2 ~ 1 - ^2 J C(,)^ U (12) 

Using Eq. 11 and 12 in Eq. 6 or 7, then allows for calculation of the echo amplitude averaged over a narrow 
energy band. A comparison between a semiclassical calculation of the echo amplitude using Eq. 11 and 12 to the full 
calculation of the previous section for the shift perturbation of a hard-wall Bunimovich billiard, is shown in fig 7. It 
is seen that the semiclassical result does not resemble the echo of an individual eigenstate, but shows good agreement 
with an average of 15 adjacent eigenstates. Since "hot" atoms are considered the population of cigenstates is a slowly 
varying function of the eigenstate number, thereby enabling the ensemble-averaged echo signal to be computed. The 
evaluation of C (io) for the hard wall billiard is done following Refs. [10] [20] [21]. The revivals occur at times ~ 2tm 
with tm being the ballistic time. 

The shift conserves the eigenspectrum of the billiard, hence mixing to nearby states is suppressed [9]. For a more 
generic perturbation such mixing effects can occur, thereby drastically reducing the range of perturbation where the 
approximation (m-f | m^} ~ 1 — is valid [10] [9]. However, the major part of the probability will still be contained 
in a very narrow "core" and for short timescales this core will not have time to dephasc. Therefore we can neglect 



where 



the decay of the core and replace |(mj | m{)f in Eq. 6 with ^| | m± 
states outside the core, and still get a reasonable estimation of the echo for short times. 



tp+ ) is the core state, sum only over 



VI. THE LIGHT-SHEET WEDGE BILLIARD 



Now the attention is turned to an experimental atom-optics billiard, namely the light-sheet wedge [6] [14]. We 
consider a wedge made from two elliptical Gaussian beams with widths of 10 [ira by 150 fim. The wedge points 
downwards so gravity confines atoms in it. The wedge angle is 98 degree and the sheets cross 67 \xm from the centers. 
Only two dimensions are considered, so the ballistic time in the third dimension is assumed to be much longer than 
the experiment time (as in Ref. [4]). Alternatively, the atoms can be confined in the third dimension by a very far 
off resonance one dimensional optical lattice, that essentially freezes the motion in that direction [7]. A long classical 
trajectory is computed for a S5 Rb atom with an energy of 1.1 x 10~ 28 J thus reaching a maximal height of 80 /im above 
the bottom of the wedge. The laser power of the wedge is taken to be 10 mW and the trap laser wavelength 779 nm 
(the resonant wavelength of the Z^-line is 780 nm). In fig 8 a Poincare surface of section is presented showing that 
the classical motion of the atom in this billiard is indeed chaotic in spite of it not being the idealized wedge billiard 
[22]. 

The two states considered \F = 2, mp = 0) (|J.)) and \F — 3, mp = 0) (||)) are separated by the hyperfine splitting 
of 3 GHz. Due to the difference in strength of dipole potential iff and ifj are related iff = fif + elldipoie with 
e = 6 x 10~ 3 for the parameters considered. 

If the wedge had been made from exponentially decaying evanescent waves close to a prism surface instead of 
Gaussian light-sheets then the change in strength in dipole potential can be described as a shift of the entire potential 
along the vertical axis. Since the shift perturbation conserves the eigenspectrum such traps are expected to yield 
very long coherence times and the use of echo techniques would not be necessary [9]. However, for a trap made 
from Gaussian light-sheets as we here consider, the perturbation does not have a simple geometric interpretation 
and does not conserve the eigenspectrum. Therefore echo techniques is expected to be of relevance, and the effects 
demonstrated is not expected to be limited to eigenspectrum conserving perturbations. 

Using the classical trajectory the echo amplitude is predicted by the methods of the previous section. The results 
are shown in Fig. 9 (|(m-f | m^)| 6 ~ 0.96). In analogy to the hard-wall billiard example given above, the echo is 
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seen to revive at times around the ballistic time. This is the same phenomena as was observed in [4] for atoms in a 
Gaussian trap. These results indicates that the revivals demonstrated in [4], are not due to the Gaussian trap being 
nearly harmonic, but a more generic feature of dipole traps even when the underlying dynamics is chaotic. 



VII. CONCLUSIONS 



A perturbative expression for the micro wave echo coherence for optically trapped atoms was derived, and confirmed 
numerically, showing that the echo signal can be written as a function of the survival probability or fidelity of 
eigenstates of the motion of the atoms in the trap. Based on this perturbative result the long time echo coherence of 
atomic ensembles was predicted. It was demonstrated that in the perturbative limit the echo signal and the survival 
probability display "system specific" features, despite the fact that the underlying classical dynamics is chaotic. 
This was done numerically by demonstrating partial revivals in the echo signal and the survival probability for a 
perturbation consisting of a small shift of the potential. Based on the perturbative expression for the echo signal 
a semi-classical method for calculating the averaged echo signal was presented. The validity of this expression was 
confirmed numerically for a hard wall Burnimovich billiard. The method was then used to calculate the echo signal 
for atoms in an experimentally realizable atom-optics billiard, namely the light-sheet wedge billiard. Revivals in the 
echo coherence were found to appear in this system, leading to the conclusion that the revivals demonstrated in [4], 
are not limited to nearly harmonic traps, but are a generic feature of dipole traps. 
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FIG. 1. The matrix elements for MW-transitions depend on the overlap of the external motion states. The trap shown is a 
Gaussian trap, and gravity gives the slope. 
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FIG. 2. Highly excited eigenstate in a Burnimovich billiard (excitation number ~ 11000). 
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FIG. 3. I (fc| I ri|) 1 2 (the elements of the transform matrix T squared (see text)) for 700 states centered around n = m ~ 11000 
for the Bunimovich stadium of Fig. 2, perturbed by a shift of 0.004. Distinct sidebands are observed. 
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FIG. 4. a: Echo amplitude from two eigenstates with similar energies in a Burnimovich billiard (excitation number ~ 11000), 
and the perturbation being a shift of 0.008. Clear revivals are seen, b: The same behavior is seen for the survival probability 
as expected from Eq. 7. 
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5. a: Solid line: Echo amplitude calculated by full calculation for a shift of 0.001 of a Bunimovich stadium, x: 
approximation by Eq. 7. Excellent agreement is seen, b: The same as (a) but for a shift of 0.008. Good agreement still exist 
in spite of very strong oscillations, c: RMS deviation between exact Echo amplitude and perturbative approximation for time 
0-0.1 as a function of perturbation (x). It is seen to agree with e 3 behavior (solid line). 
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FIG. 6. The difference between echo amplitude and its perturbative approximation as a function of time for a shift of 0.001 
of the Bunimovich Stadium. No divergence is observed even for long times. 
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FIG. 7. Dash-dotted line: Echo amplitude of single eigenstate for a Bunimovich stadium and a shift of 0.001. Dashed line: 
Semiclassical approximation using Eq. 11 and 12 in Eq. 6. Solid line: Ensemble average of echo amplitudes from 15 adjacent 
eigenstates. The semiclassical approximation does not resemble the echo from a single eigenstate, but closely resembles the 
average of just 15 eigenstates. No fitting or scaling is used. 
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FIG. 8. Poincare surface of section for the light-sheet wedge billiard described in the text. The vertical position is plotted 
as a function of the vertical momentum, when the atom crosses the symmetry axis of the wedge, for a single long trajectory. 
Phase space is seen to be chaotic. 
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FIG. 9. Echo coherence for a light sheet wedge billiard described in the text. In spite of the underlying classical dynamics 
being chaotic oscillations are seen. 
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